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Optically trapped nanoparticles have recently emerged as exciting candidates for tests of quan¬ 
tum mechanics at the macroscale and as versatile platforms for ultrasensitive metrology. Recent 
experiments have demonstrated parametric feedback cooling, nonequilibrium physics, and temper¬ 
ature detection, all in the classical regime. Here we provide the first quantum model for trapped 
nanoparticle cooling and force sensing. In contrast to existing theories, our work indicates that 
the nanomechanical ground state may be prepared without using an optical resonator; that the 
cooling mechanism corresponds to nonlinear friction; and that the energy loss during cooling is non¬ 
exponential in time. Our results show excellent agreement with experimental data in the classical 
limit, and constitute an underlying theoretical framework for experiments aiming at ground state 
preparation. Our theory also addresses the optimization of, and the fundamental quantum limit to, 
force sensing, thus providing theoretical direction to ongoing searches for ultra-weak forces using 


levitated nanoparticles. 

PACS numbers: 42.50.-p, 42.50.Wk, 37.10.Vz, 62.25.- 

INTRODUCTION 

Optically trapped nanoparticles can support explo¬ 
rations of macroscopic quantum mechanics as well as ul¬ 
trasensitive metrology very well since they can be isolated 
from the environment in a trap, cooled, and detected - 
all using a single laser beam without the need for an op¬ 
tical cavity Experiments with optically trapped 

harmonically oscillating subwavelength dielectric parti¬ 
cles IBHTU] (see Fig. have recently realized feedback 
cooling mm. nonlinear dynamics [5], non-equilibrium 
physics jl], coupling to spin degrees of freedom m and 
thermometry [3] . All experiments thus far have been car¬ 
ried out in the classical regime. 

While several groups are currently exploring ways to 
access the nonclassical regime of such systems, we present 
here the first quantum theory of trapped nanoparticle op¬ 
tical feedback cooling and force sensing. The impetus for 
investigating cooling comes from the fact that although 
levitated particles have been successfully cooled in optical 
resonators (to 10 K [T3] and 64 K M). cavityless cooling 
has been able to reach much lower temperatures (50 mK 
m) starting from the same initial (room) temperature. 

It is however an open theoretical question as to whether 
the quantum ground state can be prepared without using 
a cavity. Our model addresses this important question 
and reveals a viable route to nanoparticle ground state 
preparation for ongoing cavityless experiments. It also 
identifies the dissipative mechanism underlying experi- 


g 

mental cooling to be nonlinear in nature, in contrast to 
standard experimental techniques, which depend on lin¬ 
ear damping, and standard theory, which relies on linear 
response analysis Lastly, the model yields an¬ 

alytical results for the ensuing nonexponential decay of 
phonon number, which shows excellent agreement with 
our experimental data in the classical regime. The moti¬ 
vation for investigating force detection is provided by the 
use of levitated cavityless nanoparticles in several ongo¬ 
ing searches for various ultraweak forces [3 . Again 

it is an open question as to whether and to what extent 
force detection is limited by the effects of quantum back- 
action in those systems. To clarify this issue we derive 
in this article expressions for the force sensitivity and 
the standard quantum limit of force detection. We ex¬ 
pect our new results on cooling and force sensing in the 
quantum limit will be invaluable as this still nascent field 
matures. 

A schematic of the physical system under considera¬ 
tion is shown in Fig. [l] A subwavelength polarizable 
dielectric sphere is confined at the focus of a Gaussian 
trapping beam, and its motion is detected using a probe 
beam, polarized orthogonal to the trap. The detected 
signal is processed and fed back to the trap beam to cool 
the particle. We analyze this configuration by dividing 
it into a ‘system’ and a ‘bath’. The system consists of 
the nanomechanical oscillator and the optical probe and 
trap. The bath consists of the optical modes into which 
light is scattered by the nanosphere, and the background 



FIG. 1: a) Image of the trapped nanoparticle, b) Schematic of the experiment modeled in this article. 


thermal gas present in the experiment. We proceed to 
derive a quantum model by identifying the electromag¬ 
netic modes relevant to the problem, constructing the 
system and bath Hamiltonians, and deriving the master 
equation for the system [24l [23. All conclusions stated 
in this article follow from this master equation. 

MODEL 

The configuration Hamiltonian can be written as, 

H = Hm + Hf + Hint. ( 1 ) 

In Eq. Q, the first term on the right hand side represents 
the mechanical kinetic energy Hm — |p|^/2m, where p 
is the three dimensional momentum of the nanoparticle, 
and m its mass. The second term in Eq. 0 is the field 
energy Hf = eof |E(r)|^d^r, where E(r) is the sum of 
the trap Et, probe Ep, and background Eb electric fields. 
We model the trap and probe modes as Gaussian beams, 
and the background using a plane wave expansion. We 
find after dropping a constant term, 

Hf = hajpO^ a + E d^k;ii:jka]^(k)a;x(k), (2) 

which is simply the sum of the probe and background 
field energies with a and a^(k) representing the corre¬ 
sponding standard bosonic field operators. Finally, in 
Eq. 0 the interaction Hamiltonian is given by iJint = 
— Jy P(r) • E(r) d^r/2 = — Op Jy |E(r)|^ d^r/2, where we 
have assumed that the dielectric has a volume V, and 
that it has a linear polarizability density Op, i.e. the 
polarization density is P(r) = apE(r). Using the ex¬ 
pressions for the electric fields from the Supplementary 
Material we can evaluate iJint for small particle displace¬ 
ments r, and rewrite Eq. 0 as 

H = Hs + Hb+Hsb, (3) 


where this system Hamiltonian is 

Hs = hujpa'^a + ^ fvuJjh^^bj — ^ hgja)a(bj -!-&]), ( 4 ) 
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with mechanical trapping frequencies ujj , optomechanical 
coupling constants gj, and mechanical operators which 
obey the standard commutation relations = 1 

(j = {x,y,z}). In Eq. 0 , the bath Hamiltonian is 
Hb = d^k/ia;kaj^(k)ap(k) and the system-bath in¬ 

teraction Hamiltonian is Hsb = —CcCo /yd^r[Et(r) -|- 
Ep(r)]-Eb(r), which represents the scattering of the trap 
and probe fields into the background (see Supplementary 
Material for details). 

We now trace over the bath modes, applying the stan¬ 
dard Born and Markov approximations, since the system- 
bath coupling is weak and the bath correlations decay 
quickly [mill]. We also trace over the x and y degrees of 
particle motion, since the dynamics along the three axes 
are independent of each other, and it suffices to analyze a 
single direction m- The net result of our calculation is 
a master equation for the density matrix p{t) describing 
the optical probe and the ^-motion of the nanoparticle 

P(i) = - YT>[Q]p-b/:scP, (5) 

where the first term on the right hand side represents 
unitary evolution of the system with Hg = hujpO^a -\- 
hu3zb\bz — hgza^a{bz + bl). The second term corresponds 
to the positional decoherence of the nanoparticle due 
to scattering of trap photons, with the Lindblad super¬ 
operator 'D[Qz]p = QIQzP + pQIQz - 2QzpQl, where 
Qz = b\ + bz, and At is the heating rate due to trap beam 
scattering as defined in the Supplementary Material. The 
third superoperator describes the loss of photons from 
the probe, also due to scattering by the nanoparticle, 
'^sc[p(t)] = -B (v[a\ + {7uj'^el/5c^)T>[aQz]^ p, where 4 
is the oscillator length. 
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The nanoparticle also experiences collisions with back¬ 
ground gas particles at the ambient temperature T. This 
effect may be accounted for by adding to the right hand 
of Eq. ([^ the superoperator [26] 

B[p{t)]=-^V[Q,]p-^V[P^]p 

2 2 (6) 

“ ’ 

where — bz), and curly braces denote an anti¬ 

commutator. The first term on the right hand side corre¬ 
sponds to momentum diffusion and Dp = 2pfkBT£^/h^, 
where ks is Boltzmann’s constant. The second term de¬ 
scribes position diffusion with Dq = pfh?/{2AkBTm?£‘l). 
The third term accounts for friction, and by Stokes law 
we have rjf = Grr/rr^, where is the radius of the 
nanoparticle and p is the dynamic viscosity of the back¬ 
ground gas. As shown earlier, internal and center-of- 
mass heating of the nanoparticle due to optical absorp¬ 
tion and blackbody radiation are negligible in systems 
such as ours, as are particle size and shape effects, as 
well as trap beam shot noise [IIIEt]. 

We now characterize the measurement of the oscillator 
displacement using input-output theory from quantum 
optics |55| applied to the nanoparticle. Specifically, the 
incoming probe field Uin interacts with the nanoparticle, 
and the outgoing probe field Oout carries a signature of 
this interaction (as shown in the Supplementary Mate¬ 
rial) 

^out — ^in T (7) 

where y = ApzAt is the scaled optomechanical coupling, 
with integration time At (determined by the detection 
bandwidth), and we have written the probe beam as a 
coherent state a = —ia + v, with a a classical number 
and V a bosonic annihilation operator. A homodyne mea¬ 
surement on the output field yields a current [3S] 



FIG. 2: a) Shows the diffusive evolution of the trapped 
nanoparticle’s position at atmospheric pressure, b) Har¬ 
monic motion of the nanoparticle at a lower pressure of 
4 X 10“^ mbar. The reduction of the amplitude of harmonic 
motion corresponds to the turn-on of feedback, i.e. cooling. 

and the feedback current is /fb = {Pz) + {t), 

where ^'(t) has the same properties as ^{t). This 
form of the Hamiltonian implies a feedback force Ffb = 
—dHfi^/dQz which is equivalent to that used in experi¬ 
ments in the classical regime [Sj. Taking the Markovian 
limit where the feedback occurs faster than any system 
timescale, and applying quantum feedback theory for ho¬ 
modyne detection [30] . we find that the following super¬ 
operator must be added to Eq. (§ 


h = x^m.){t) + v^m, ( 8 ) 


where <P = a^Aut is the average detected flux of probe 
photons, and ^(t) is a stochastic variable with mean 
= 0 and correlation {^{t)^{t')) = S(t — t'). 

In the experiment, the detected current Ih is frequency 
doubled, phase shifted, and fed back to modulate the 
power of the trapping beam m- This results in a feed¬ 
back Hamiltonian i/fb = where G is the di¬ 

mensionless feedback gain related to the trap intensity 
modulation [511^. 


It OJz 


(9) 


J^[p{t)] = -tx^<l>G[Ql{Pz,p}] - ^G^V[Ql]p, (10) 

where the first term on the right hand side represents 
the desired cooling effect of the feedback, and the second 
term the accompanying backaction. We emphasize that 
in contrast to standard optomechanics, the feedback and 
backaction terms are highly nonlinear in the oscillator 
variables. The presence of this nonlinearity distinguishes 
our system from conventional cavity optomechanics and 
results in qualitatively different dynamics, as we show 
below. 

The full master equation, assembled from Eqs. (§, (§ 
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FIG. 3: a) The y and z phonon cooling dynamics [Eq. (13l]. b) Steady state phonon number versus pressure [Eq. (14al]. 
Circles represent experimental data, and the solid curve our theoretical model for a fused silica sphere (cr = 2.1 and density = 
2200kg/m^) radius Vd = 50nm, 1064nm trap (lOOmW) and probe (lOmW) beams, a mechanical frequency a;z/27r = 38kHz, 
y Rs 10“^, and trap intensity modulation M < 0.1%. The dotted lines represent the equivalent curves for one of the transverse 
degrees of freedom {ijjy/2'K = 138 kHz). The dashed curve in a) represents the prediction of our theory for a setup placed in a 
cryostat with the feedback chosen optimally, keeping M < 10%. 


and (101 is then 


m = - (A + ^p)I?[Q,]p(t)/2 

- BV[a]p{t) + B[p{t)] + T[p{t)], 


( 11 ) 


where the new system Hamiltonian Hs = hujpV^v + 
JTij->zb\bz — iahgziv^ — v){hz + bl), accounts for the lin¬ 
earization of the probe implemented above, and the mas¬ 
ter equation now includes a mechanical decoherence term 
due to scattering from the probe, in addition to the trap 
beam, with heating rate Ap. 

A sample experimental data set of the measured po¬ 
sition of the nanoparticle along the y axis is shown in 
Fig.i At atmospheric pressures. Brownian effects, as 
given by Eq. ([^ dominate, and the position of the parti¬ 
cle follows a diffusive evolution, as can be seen in Fig.[^. 
At lower pressures the particle’s evolution becomes in¬ 
creasingly ballistic. The ensuing harmonic motion is 
shown in Fig. @3, both in the absence as well as in pres¬ 
ence of feedback [Eq. @]- The decrease in amplitude of 
the harmonic motion is due to the presence of parametric 
feedback cooling. 


PHONON DYNAMICS 

Employing the master equation [Eq. ( 0 ] to consider 
the question of ground state occupation, tracing out the 
optical probe field, and using the resulting reduced mas¬ 
ter equation for the nanoparticle only, we find the equa¬ 
tion for the dynamics of the phonon number {N = blbz), 
(N) = - J(N^) - K{N)+L, where J =[12G-5AG^]X^$, 


K = pf /2m + J, L = D — J/2, the dot denotes a time 
derivative, and D = + Dg with D/ = Dp + At + Ap 

accounting for positional decoherence. We assume that 
the nanoparticle is described by a thermal state EUSl- 
IM] . for which (A^) = 2{N)'^ + (N) [^, a relation which 
simplifies the phonon dynamics to 

{N)=-2J{Nf-{J + K){N)+L. (12) 


The total effect of parametric feedback on the phonon 
dynamics is contained in the parameter J, which is de¬ 
termined by the difference between the feedback cooling 
and backaction heating. In the experiments J 7 ^ omm, 
making the phonon dynamics of cooling nonlinear, and 
the oscillator energy loss nonexponential, as shown be¬ 
low. We stress that this behavior is qualitatively dif¬ 
ferent from standard quantum cavity optomechanical 
theory, which characterizes cooling as a linear damp¬ 
ing process resulting in an exponential decay of energy 
(see, e.g. [m [32] and Eq. (82) in [H]). We note that 
G = Gopt = 1/9 maximizes (N) in Eq. (12) with the 
maximum nonlinear cooling rate Jmax = 2x^^/3. 

Assuming the initial condition {N{0)) = Nq = 

D'p2m/pf = kBTeff/huJz, where T^g is the effective tem¬ 
perature of a bath due to gas and optical scattering com¬ 


bined, the analytical solution to Eq. (12) is 


(Nit)) = - 


(J + K) 
4J 


-tanh I —hi 

2Jt \t 


(13) 


where 9 = tanh ^ [{2JNo + J + Ar)r] 
timescale r = 2 [(J -h AT)^ -h 8JL] 


and the cooling 


From Eq. (13) 
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the steady state phonon number is 


10 


N.. s ^Ito (IV(i)) = ^ 


Vf No 
2m 2 J 


Dp -\- At -\- Ap 


2J 


(14a) 

(14b) 


where the approximation is valid for Nq ;g> 1. To reach 
the ground state, we need to maximize the feedback cool¬ 
ing J, which can be done by setting G = Gopt- We also 
need to minimize gas heating, which can be accomplished 
by going to low pressures and cryogenic temperatures, 


such that Dp is negligible in Eq. (14). Below we discuss 


situations involving realistic experimental parameters. 

Two plots of the nonlinear phonon dynamics are shown 
in Fig. ir for the z and y motion at 10 ^ mbar along 
with experimental values measured by us (circles). The 
solid curve represents {N{t)) as given in Eq. (13), while 
the dotted curve gives the corresponding equation for 
motion along one of the transverse directions y (which is 
nearly degenerate with x i.e. ujy « Wx), see Fig. [^. In 
Fig. Hr we show three plots of the steady state phonon 
number, as the vacuum pressure is tuned. The solid and 


dotted curves represent Nss [Eq. (14a)] at 300 K for z and 
y motion respectively, while the circles are experimental 
data. As can be seen, in all cases there is very good 
agreement between theory and experiment. The dashed 
curve in Fig. predicts the steady state phonon number 
for an identical configuration, but placed in a cryostat at 
4K. The ground state can be prepared if starting at 
high pressures, the particle is cooled while continuously 
increasing the feedback gain as in in, and keeping the 
trap modulation M = 10%. Proceeding in this manner, 
we find that below < 10“® mbar optimal feedback J = 
Anax can be achieved, and the ground state occupied. 

We note that practical cooling to lower phonon num¬ 
bers is currently limited by a number of factors. These in¬ 
clude high pressures enforced by nanosphere loading tech¬ 
nologies, classical errors from the electronic feedback loop 
and laser noise, measurement uncertainties due to detec¬ 
tor bandwidth limitations, and collection inefficiencies of 
the scattered light [10 El [H]. However, these prob¬ 
lems are technical rather than fundamental, and efforts 
are underway to overcome these limitations [S5] , suggest¬ 
ing strongly the possibility of cooling to the ground state 
using the parameters presented in this article. 


FORCE SENSING 

We now consider force sensing using the nanoparticle 
model given by the master equation of Eq. 0- Since 
the state of the nanoparticle is continuously monitored, 



uj /2 t : (kHz) 


FIG. 4: Experimentally measured positional PSDs for all 
three degrees of freedom with dark lines representing the 
theoretical fits to the data [Eq. ( |18[ )]. Data was taken at a 
moderate vacuum pressure of 10 mbar and clearly shows the 
Lorentzian shape of the resonance. These fits were used to 
extract the values of u)j, F, and Nsa- 


the master equation can be unraveled in terms of a set 
of Langevin equations describing the evolution of the 
quadratures Qz and Pz plus a stochastic force due to 
the measurement backaction which gives [371 [3H] 

Qz = C*[Qz]=ujzP. 

Pz = Ct[Pz]+F/miOziz (15) 

= -tOzQz - PPz + F/muJziz, 

where is the Liouvillian superoperator (dual to the 
superoperator Cq appearing in the master equation p = 
Cq[p]) defined by Tr(p£^[A]) = Tr(A£o[p]) for any arbi¬ 
trary operator A [39] . The parameter F = Ig-I-JT', where 
Fq is the gas damping and SF « 12x‘^^G{{N) -1-1/2) is 
the nonlinear feedback damping m O E]. Finally F = 
Ft + Fp is the sum of the (independent) stochastic forces 
due to thermal and feedback backaction heating, respec¬ 
tively, with zero mean and correlations {FT{t)FT{t')) = 
SpSit — t') and (FF{t)FF{f)) = SpSit — t'), with 


St = 2TO/ofc_B7eff 


Sf = 27mhujzX^^G'^ (2 {Nf + 2 (N) + l) . 


(16) 


The presence of the (Al)-dependent factor in Eq. (16) im¬ 


plies that the feedback noise is dependent on the system 
state, and is therefore non-additive. Furthermore, the de¬ 
pendence is nonlinear in {N). Both of these features are 
fundamentally different from the typical additive feed- 
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back noise in standard cavity optomechanics, which is 
independent of the state of the system [55] , 


We convert Eq. (15) into the second order differential 
equation for the position gz + = F/m, and 

take its Fourier transform to find the position spectrum 
= Xm(w)F(a;), where 


Xm(w) = < m 




— iuiF 


(17) 


is the optomechanical susceptibility of our oscillator. Fi¬ 
nally the positional power spectral density (PSD) noise 
spectrum is given by 


92(w)|^) =IXm|^ (St+Sf) + 


x2<?’ 


(18) 


where the last term in the equation comes from the shot 
noise of the measured signal [Eq. ([^]. A typical example 
data set of the positional PSDs at moderate vacuum is 
shown in Fi g.[^ along with fits to the theoretical expres¬ 
sion of Eq. (flsF. 


In view of the fact that trapped nanoparticles offer 
the possibility of ultrasensitive force measurements [51 
I21H23| , we express our measurement noise spectrum 
[Eq. @[ in terms of the estimator q(u))/xm in order 
to investigate the fundamental limits of such measure¬ 
ments. The sensitivity of force estimation is set by the 
force noise PSD 


E(w)| )= St + Sf+ Ss{u:), 


(19) 


where Ss{oj) = Ss{0) (l — (ui/ojz)'^) + {uiFjojl)'^ and 

'S's(O) = {m£z‘jjl)‘^/x^^- Only the last term carries an 
w dependence in Eq. (19). A plot of Ss{uj) is shown in 
Figj^in the high as well as low total damping F regimes, 
both of which are experimentally accessible m ng. 
The minimum value of Ss(uj), and therefore the opti¬ 
mal force sensitivity, occurs at the response frequency 
ccopt = V^^z - OV2. 


The first two terms in Eq. (19) scale linearly with 


the optical power while the shot noise scales inversely 
(i.e. St + Sf ~ ^ and Ss ^ 1/^)- Therefore there 
is a power that minimizes the total noise, representing 
the standard quantum limit for our system. Assuming 
that the feedback is optimal (i.e. J = Jmax), the stan¬ 
dard quantum limit is reached when x^^SQL ~ F/8 and 
equals 


i^r 


SQL 


2mFQkBT + AmhuJz (^At + V^F^j . ( 20 ) 


The first term in (20) represents a thermal contribu¬ 


tion from the background gas; the second term is due 



FIG. 5: a) Plots of the shot noise force PSD [the last term of 
Eq. (19l] versus the normalized mechanical frequency lu/ujz 
for low and high total damping F. The minimum occurs 
for uJopt = — F'^/2. b) Plot of the force sensitivity as 

a function of the normalized optical power ^/^sql at high 
vacuum. The standard quantum limit is reached whe n th e 
shot noise balances the recoil and backaction noises [Eq. (19)[. 


to scattering of photons from the trapping beam; and 
the third term contains the effects of light scattering 
and shot noise from the probe, as well as the feedback 
backaction. At the low vacuum pressures currently avail¬ 
able (i.e. < 10“^ mbar) the gas contribution is negligible, 

implying a minimum force sensitivity of \J (|E| )gQL ~ 
10 “^^N/-\/Hz = zN/-\/Hz and optimal probe power of 
huip<l>sQL ~ 2mW, where the remaining system parame¬ 
ters have been taken from the caption of Fig. Even at 
this limit the system can be readily used to test for vio¬ 
lations of Newtonian gravity (~ 10“^® N) [^ with mod¬ 
erate measurement bandwidths. However, backaction 
effects will impose long interrogation times on experi¬ 
ments searching for new small scale (~ 10“^^ N) [21], and 
Casimir forces (~ 10“^^N) [51 f^. Conversely for short 
measurement times, our calculations show that backac¬ 
tion effects, which are of interest in their own right in 
optomechanics m, can be observed at moderate laser 
powers and readily attainable vacuum pressures. 


6 






















CONCLUSIONS 

To conclude, we have presented a quantum model that 
describes the cooling and force sensing characteristics of 
an optically trapped subwavelength dielectric particle. 
We have shown that the predictions of this model for 
cooling are in very good agreement with experimentally 
measured occupation values in the classical regime. Fur¬ 
ther, we have demonstrated that quantum ground state 
preparation is challenging, but achievable in anticipated 
experiments. Finally, we have derived the standard quan¬ 
tum limit to force sensing, indicating experiments where 
the role of quantum backaction needs to be accounted 
for. The model presented by us opens the door to the 
characterization of the quantum behavior of a system im¬ 
portant for macroscopic quantum mechanics mm , op¬ 
tical tweezing [2], ultrasensitive metrology [3], and non¬ 
equilibrium physics [1]. With the proper identifications, 
our theory is also applicable to electromechanical systems 
with parametric feedback m- 


ACKNOWLEDGEMENTS 

We are grateful to C. Stroud, A. Aiello, B. Zwickl, 
and S. Agarwal for useful discussions. This material is 
based upon work supported by the Office of Naval Re¬ 
search under Award Nos. N00014-14-1-0803 and N00014- 
14-1-0442. ANV thanks the Institute of Optics for sup¬ 
port. LPN acknowledges support from a University of 
Rochester Messersmith fellowship. 


* Electronic address: brandon.rodenburg@gmail.com 

[1] Neukirch, L. P., Gieseler, J., Quidant, R., Novotny, L. 
& Nick Vamivakas, A. Observation of nitrogen vacancy 
photoluminescence from an optically levitated nanodia¬ 
mond. Opt. Lett. 38, 2976-9 (2013). 

[2] Gieseler, J., Spasenovic, M., Novotny, L. & Quidant, 
R. Nonlinear Mode Coupling and Synchronization of a 
Vacuum-Trapped Nanoparticle. Phys. Rev. Lett. 112 , 
103603 (2014). 

[3] Millen, J., Deesuwan, T., Barker, P. & An¬ 
ders, J. Nanoscale temperature measurements us¬ 
ing non-equilibrium Brownian dynamics of a levitated 
nanosphere. Nat. Nanotechnol. 9, 425-9 (2014). 

[4] Gieseler, J., Novotny, L., Moritz, C. & Dellago, C. Non¬ 
equilibrium steady state of a driven levitated particle 
with feedback cooling. New J. Phys. IT, 045011 (2015). 

[5] Neukirch, L. P. & Vamivakas, A. N. Nano-optomechanics 
with optically levitated nanoparticles. Contemporary 
Physics 56 , 48-62 (2015). 


[6] Li, T., Kheifets, S. & Raizen, M. G. Millikelvin cooling 
of an optically trapped microsphere in vacuum. Nature 
Phys. 7 , 527-530 (2011). 

[7] Yin, Z. Q., Geraci, A. A. & Li, T. C. Optomechanics of 
levitated dielectric particles. Int. J. Mod. Phys. B 27 , 
1330018 (2013). 

[8] Bateman, J., Nimmrichter, S., Hornberger, K. & Ul- 
bricht, H. Near-held interferometry of a free-falling 
nanoparticle from a point-like source. Nat. Commun. 5 
(2013). 

[9] Arita, Y., Mazilu, M. & Dholakia, K. Laser-induced rota¬ 
tion and cooling of a trapped microgyroscope in vacuum. 
Nat. Commun. 4, 2374 (2013). 

[10] Scala, M., Kim, M. S., Morley, G. W., Barker, P. F. 
& Bose, S. Matter-wave interferometry of a levitated 
thermal nano-oscillator induced and probed by a spin. 
Phys. Rev. Lett. Ill, 180403 (2013). 

[11] Gieseler, J., Deutsch, B., Quidant, R. & Novotny, L. Sub- 
Kelvin Parametric Feedback Cooling of a Laser-Trapped 
Nanoparticle. Phys. Rev. Lett. 109 , 103603 (2012). 

[12] Neukirch, L. P., von Haartman, E., Rosenholm, J. M. & 
Nick Vamivakas, A. Multi-dimensional single-spin nano- 
optomechanics with a levitated nanodiamond. Nat. Pho¬ 
tonics (2015). 

[13] Millen, J., Fonseca, P. Z. G., Mavrogordatos, T., Mon- 
teiro, T. S. & Barker, P. F. Cavity Cooling a Single 
Charged Levitated Nanosphere. Phys. Rev. Lett. 114 , 
123602 (2015). 

[14] Kiesel, N. et al. Cavity cooling of an optically levitated 
submicron particle. PNAS 110 , 14180-5 (2013). 

[15] Kippenberg, T. J. & Vahala, K. J. Cavity opto¬ 
mechanics. Opt. Exp. 15 , 17172-17205 (2007). 

[16] Marquardt, F. & Girvin, S. M. Optomechanics (a brief 
review). Physics 2 , 40 (2009). 

[17] Romero-lsart, O. et al. Large Quantum Superpositions 
and Interference of Massive Nanometer-Sized Objects. 
Phys. Rev. Lett. 107 , 020405 (2011). 

[18] Aspelmeyer, M., Kippenberg, T. J. & Marquardt, F. 
Cavity optomechanics. Rev. Mod. Phys. 86, 1391-1452 
(2014). 

[19] Meystre, P. A short walk through quantum optomechan¬ 
ics. Annalen der Physik 525 , 215-233 (2013). 

[20] Asenbaum, P., Kuhn, S., Nimmrichter, S., Sezer, U. & 
Arndt, M. Cavity cooling of free silicon nanoparticles in 
high vacuum. Nat. Nanotechnol. 4 , 2743 (2013). 

[21] Geraci, A. A., Papp, S. B. & Kitching, J. Short-Range 
Force Detection Using Optically Cooled Levitated Micro¬ 
spheres. Phys. Rev. Lett. 105 , 101101 (2010). 

[22] Ranjit, G., Atherton, D. P., Stutz, J. H., Cunningham, 
M. & Geraci, A. A. Attonewton force detection using 
microspheres in a dual-beam optical trap in high vacuum. 
Phys. Rev. A 91 , 051805 (2015). 

[23] Moore, D. C., Rider, A. D. & Gratta, G. Search for 
Millicharged Particles Using Optically Levitated Micro¬ 
spheres. Phys. Rev. Lett. 113 , 251801 (2014). 

[24] Pflanzer, A. C., Romero-lsart, O. & Cirac, J. I. Master- 
equation approach to optomechanics with arbitrary di¬ 
electrics. Phys. Rev. A 86, 013802 (2012). 

[25] Carmichael, H. J. Statistical Methods in Quantum Op- 


7 



tics 1: Master Equations and Fokker-Planck Equations 
(Springer, 2002). 

[26] Diosi, L. Quantum Master Equation of a Particle in a Gas 
Environment. EuroPhys. Lett. (EPL) 30, 63-68 (1995). 

[27] Chang, D. E. et al. Cavity opto-mechanics using an opti¬ 
cally levitated nanosphere. PNAS 107, 1005-10 (2010). 

[28] Gardiner, C. W. & Zoller, P. Quantum Noise: A 
Handbook of Markovian and Non-Markovian Quantum 
Stochastic Methods with Applications to Quantum Optics 
(Springer, 2004), 3rd edn. 

[29] Mancini, S., Vitali, D. & Tombesi, P. Optomechanical 
Cooling of a Macroscopic Oscillator by Homodyne Feed¬ 
back. Phys. Rev. Lett. 80, 688-691 (1998). 

[30] Wiseman, H. M. & Milburn, G. J. Quantum theory of op¬ 
tical feedback via homodyne detection. Phys. Rev. Lett. 
70, 548 (1993). 

[31] Wilson-Rae, I., Nooshi, N., Zwerger, W. & Kippenberg, 
T. J. Theory of Ground State Cooling of a Mechanical 
Oscillator Using Dynamical Backaction. Phys. Rev. Lett. 
99, 093901 (2007). 

[32] Marquardt, F., Chen, J. P., Clerk, A. A. & Girvin, S. M. 
Quantum Theory of Cavity-Assisted Sideband Cooling of 
Mechanical Motion. Physical Review Letters 99, 093902 
(2007). 

[33] Genes, C., Vitali, D., Tombesi, P., Gigan, S. & As- 
pelmeyer, M. Ground-state cooling of a micromechanical 


oscillator: Comparing cold damping and cavity-assisted 
cooling schemes. Phys. Rev. A 77, 033804 (2008). 

[34] Romero-Isart, O. et al. Optically levitating dielectrics in 
the quantum regime: Theory and protocols. Phys. Rev. A 
83, 013803 (2011). 

[35] Gerry, C. & Knight, P. Introductory Quantum Optics 
(Cambridge University Press, Cambridge, 2004). 

[36] Mestres, P. et al. Cooling and manipulation of a levitated 
nanoparticle with an optical fiber trap. Appl. Phys. Lett. 
107, 151102 (2015). 

[37] Gisin, N. & Percival, I. C. The quantum-state diffusion 
model applied to open systems. J. Phys. A. Math. Gen. 
25, 5677-5691 (1992). 

[38] Halliwell, J. & Zoupas, A. Quantum state diffusion, den¬ 
sity matrix diagonalization, and decoherent histories: A 
model. Phys. Rev. D 52, 7294 (1995). 

[39] Hornberger, K. Entanglement and Decoherence, vol. 768 
of Lecture Notes in Physics (Springer Berlin Heidelberg, 
Berlin, Heidelberg, 2009). 

[40] Purdy, T. P., Peterson, R. W. & Regal, C. A. Observation 
of radiation pressure shot noise on a macroscopic object. 
Science 339, 801-4 (2013). 

[41] Villanueva, L. G. et al. A Nanoscale Parametric Feedback 
Oscillator. Nano Lett. 11, 5054-5059 (2011). 



Supplemental Materials: Quantum Model of Cooling and Force Sensing With an 

Optically Trapped Nanoparticle 


ELECTRIC FIELDS 


The total electric field as defined in the main article is written as 

E(r,t) = Et(r,f) + Ep(r,t) + Eb(r,<). 


(SI) 


The Et and Eb helds are both assumed to be Gaussian beams with the trap treated classically and the probe treated 
as a quantized beam with frequency ujp, linewidth Alu, waist wq, and canonical bosonic operators [a,a^] = 1 |S1) 


Ep(r,t) = i 


( hujpAuj \ 
V 47reoC ) 


e“^(^/=-*)epG(r', wp)a + h.c. 


(S2) 


where G(r, oj) is the spatial mode function of the beam. The background field Eb is simply the quantum held of all 
other modes and can be represented via a typical plane-wave expansion |S2] . 


FREE FIELD HAMILTONIAN 

Our conhguration Hamiltonian as given in the main text is 


H = Hm + Hf+Hint, (S3) 

where Hm = I2m is the particle’s kinetic energy for the momentum p = ^Jp"^ Py ^'^'4 Lfint are the 

free held and interaction Hamiltonians respectively. The energy of the free held Hamiltonian is 

Hf=eoj\E[v,t)\^d?v. (S4) 


The term proportional to |Et(r)|^ corresponds to the energy of the trap beam, and can be neglected as it represents 
a constant offset of the Hamiltonian. The energy of the free probe held is given by [ST] 

eo J IEp(r) I^ d^r = fiWpa^a. (S5) 

The energy of the background held, which can be found in several textbooks (e.g. |S2) 1 can be written in a plane wave 
expansion as 


Hb=^ d^k;iwka|,(k)a^(k). 


(S6) 


The cross term between the trap and background helds, given by J (Et -f Ep) • Ebd^r, vanishes as their mutual 
overlap is very small. This cancellation also represents the avoidance of self-interference and mode overcounting in 
our model. The cross term between the probe and the background helds vanishes for the same reason. The cross term 
between the trap and probe helds given by Et ■ Ep, vanishes due to polarization orthogonality. Finally, combining 


Eqs. (S5) and (S6), the free held Hamiltonian is 


Hf = hujpo) a + d^k;iwkajj(k)a;_.(k). 


(S7) 


which is simply the sum of the probe and background held energies. 
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INTERACTION HAMILTONIAN 


The interaction Hamiltonian between the fields and the nanoparticle is given by 

1 




P(r) • E(r) d'^r. 


(S8) 


Assuming that the particle has a linear polarizability described by P(r) = apE(r), the interaction Hamiltonian is 

EoCc 




E(r) d^r, 


(S9) 


where Cc = 3(er — l)/(er + 2) is the Clausius-Mossotti relation for the effective relative permittivity of a dielectric due 
to local field effects, and V denotes integration over the volume of the dielectric particle. 

Now when we use the total electric field from Eq. (SI I in the above equation, we again get a number of terms. The 
terms proportional to |Eb| and |Et ■ Ep| we neglect as these represent a renormalization of the background modes 
and we again assume the trap and probe are cross-polarized. 


Trap potential 


The first non-negligible term is the effect of the trap beam proportional to Jy |Ep|^ d^r. Since the nanoparticle is 
smaller than the wavelength of any relevant optical field, the integration can be written as fy — V J <5(q)d^r, 
where q is the center of mass position of the particle. Therefore iJint for the trap-particle interaction is given by 


-^H|Et(q)|" = -^HE2 


ecEo 


-^0^||G(q)r 


Qz 


VE^ I -1 + ( ^ 

ZR 


+2('^y 



Oiq% 


(SIO) 


where Qz is the longitudinal coordinate, and Qx and qy are the transverse coordinates. Ignoring the overall constant 
in the above equation, we finally have for our trap Hamiltonian 


I^trap — 


2m + S” 


P 


2 2 


9 2 

■ ‘^xqx 


2 2 
■ ^yqy 


)■ 


(Sll) 


where ojz = \Jtc.^^E'^Vj(mzj^) and u)x,y = \/j(mwQ). Writing our canonical nanoparticle variables as 


qj = a/S/ {2mujj){b^j + bj) = ^jib^A -l- bj) = ijQj and pj = iy/hmujj/2{b'^ — b) allows us to rewrite Eq. (Sll) as 


H, 


trap 


= E 


^hwjb]bj. 


Optomechanical coupling 

Now we examine the coupling between the probe beam and the trapped nanoparticle given. The Hamiltonian for 
this term is 


Hom = / |Ep(r)|^d3r 


= -V 


= -V 


V 

tchujpAui 

Anc 

CchuipAuj 

2tt‘^WqC 


+ ^ ) |G(q- Ar)| 


(S12) 


(a'^'a + + ^i^xqx + Ayqy) + C>(q^)^ , 
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where we have assumed that the probe beam is shifted from the trap beam by a small amount Ar = {Ax, Ay, Az). 
Now the terms proportional to a constant times either a^a or qj (for j G {x,y,z}) can be incorporated into shifts in 
the optical and mechanical frequencies, and the oscillator position. The optomechanical coupling term between the 
probe and particle is then given by 

Hom = ) = - E (S13) 

\ ^R J j 

The system Hamiltonian is then a combination of the energies of the probe field, optical trap, and optomechanical 
coupling and is given by 


Hs = huJpO^a + E] ~ E/ + ^])' 


(S14) 


Optical scattering 


Computing the cross coupling Hamiltonian between the field and the background due to scattering by the nanopar¬ 
ticle gives 




Et • Eb d^r 


for the trap beam and 


(1^) / E ■ ^-L + - ^o)qz) 

X (e^(k) • E5a^(k)e-*(“'=-‘"‘)‘ + h.c.) , 


-ec€o / Ep • Eb d'^r « -i 


IV 


fiecV 

Stt^wq 


J d^kyC^E ■ ^-L + “ ko)qz) 


X (ep • e^(k)aj,(k)ape h.c.^ 


(S15) 


(S16) 


for the probe beam, assuming elastic scattering and amplitude of the particle motion small relative to an optical 
wavelength [S3j . Applying standard Born-Markov theory and tracing over the background optical modes and particle 
motion transverse to the z-axis allows us to derive a master equation 


P(i) 


Sc6 eoTT 

-Y'^[Qz]+Csc[p{t)] 


1 {hz0p)^elv^ 
2 Ait^Wq{c/A u}) 



5c2 



p{t) 


(S17) 


For reference, if we do not trace over the transverse motion, the three-dimensional analogue of the trap scattering 
is 

cr [p(0] ^ 

where |j and x indicate transverse motion parallel (or not parallel) to the polarization direction of the trap beam. A 
similar expression holds for the probe 


[p{t)] 


1 {hwpyelV^ 
h'^c^ 2 Att'^Wq{c/A uj) 


2 ^2 

25[ap] + {^[aqy] + 2T>[ag||] -f TD[aqA^ 


p{t). 


(S19) 
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DERIVATION OF THE DETECTED HOMODYNE CURRENT AND FEEDBACK 


In order to find the input-output relations for the probe field, we write the system Hamiltonian (for the single 
degree of freedom Q^) as Hg = fvjjpO^a + fujjJ)\bz — Hga^aQ^. If we move into the interaction picture for the probe 
field where a —>■ then this becomes Hg = hujzb\bz — hga^aQz. 

We assume the probe field is initially a coherent state which can be written as a = —ia + v, where a is a constant 
and z; is a field annihilation operator. In this case the optomechanical coupling in our system Hamiltonian becomes 

-ffoM = —figzQzO' o, = hgzQz + iav"^ — ia*v^ . (S20) 

The term proportional to \a\^ Qz is responsible for simply shifting the mean position of the oscillator and can safely 
be ignored, leaving us with Hqm = ifigzQzictv^ — a*v). 

The Heisenberg equation of motion for v is given by 

i) = ^[a, H] = agzQzi (S21) 


which can be integrated formally to give 


v{t) = v{to) + / agzQz dt' « v{to) + agzQz^t, 
Jto 


(S22) 


where the integration is taken over a time At short compared to 1 /lOz ■ If we had picked an initial time tf > t then we 
would have computed v{t) = v(tf) — agzQzAt. By identifying the input state as Oin = —ia + v{tQ), and the output 
state as Oout = then we can relate the output and input fields by 


^out — ^in “t” ^agzAt — din T 

where we have defined the variable y = AgzAt. 

Now detected homodyne current is proportional to 

4 = x"^(0z)(t) + \/x^?W, 


(S23) 


(S24) 


where we have converted from photon numbers to rates by using the identity for the photon flux <!> = {a) a) Auj = a'^Aut, 
and have introduced the stochastic variable ^{t) which is due to the shot noise of the detection. 

As described in the main text, phase shifting the measured signal (or equivalently adding a short time delay) is 
equivalent to measuring a different quadrature of motion, i.e. the current that is fed back is /fb = (Pz) + \/x^^'(i) 
, where ^'(t) has the same properties as ^{t). Now writing cr = x^'^Pz allows us to write the master equation for the 
feedback in standard notation [S4] 


p = IC[ap + pa'’] + 


QlAPz.p} 




G^V[Ql]p 


(S25) 


where the Liouvillian superoperator K, is defined as lC[p] = [F,p\/ih and where F is the feedback term that comes 
from the feedback Hamiltonian iJfb = IfbF, and is chosen to be F = HGQ^ to match the classical physics [S5] of the 
problem as described in the main text. The gain coefficient G may be related to the experimental trap beam intensity 
modulation M = Alt/It by using the fact that Hpo = biOzMQ’/. = IfbF, and therefore [S5| 


M = 


Gx^<P{Qz){Pz) 

Wz 


Gx^^N 

UJz 


(S26) 
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FULL MASTER EQUATION 


In the main text of this Letter we presented the master equations describing the density matrix of only the single 
degree of freedom represented by the 2 motion of the nanoparticle. For convenience, we present the full master 
equation of the three-dimensional motion of the nanoparticle, as well as of the probe beam itself, which can be used 
to calculate the statistics of any quantum observable for the nanoparticle or probe beam. 

The master equation for the density matrix p describing the full three-dimensional motion of the nanoparticle, as 
well as the probe field is given by 


p{t) = — [Hs,p{t)] j Unitary dynamics 

[p(t)] + [p(i)] i Photon scattering 


(S27) 


PfksT 


q, [q,pW] 




UkeTm^ I 


p, [p,pW] 


2mh . 


q, {p,p(t)} > Gas scattering. 


This equation can be used for studying full dynamics, light-matter entanglement and photon statistics. 


HOMODYNE PHOTODETECTION: A COMPARISON OF THEORY AND EXPERIMENT 


The experimental data in Fig. 2 of the main article was obtained by using a scheme more involved, and of more 
general applicability, than indicated in the simplified detection model presented in the main article. In the laboratory, 
the position of the nanoparticle was determined from the optical interference between the unscattered probe and a 
spherical wave radiated by the induced dipole of the polarizable sphere. The interference signal varies linearly with the 
particle position for oscillation amplitudes small compared to the optical wavelength. In order to detect the particle 
position and also to eliminate the large constant background term (equal to the unscattered probe flux « (p) in an 
experimental setting, we leveraged a balanced homodyne detection scheme, as described in [S3| and |S5) . 

In this configuration, one of the detectors sampled the entire spatial profile of the probe beam, while the other was 
apodized so that only the center of the beam was sampled. As the trapped particle moved in the axial direction, 
the relative amount of scattered light collected by each detector, and thus the AC term of the resulting homodyne 
current, was modulated. The optical channels were adjusted so that each detector recorded the same average power, 
and the resulting homodyne signals were subtracted to eliminate the common DC bias. 

In our theoretical model, rather than include an additional quantized mode corresponding to the full induced dipole 
emission pattern, we considered for simplicity only the field scattered back into the probe mode by the nanoparticle. 
This simplification is admissible as long as the value of the linear optomechanical coupling constant x in th® model is 
taken from experiment |S3) . We emphasize that our experimental method of feedback cooling works even if the trap 
and probe foci are very close to each other, in which case the optomechanical coupling is essentially quadratic in the 
particle position, while the position detection (using the dipole wave mode) is still linear. In the theoretical calculation 
of the phonon number, the linear (or quadratic) optomechanical coupling simply adds a position offset (or frequency 
shift) to the trap. Both the offset as well as the shift are negligible for weak probe light used in the experiment, 
rendering the theoretical predictions of ground state cooling identical for linear as well as quadratic optomechanical 
coupling. 


COMPUTING THE GAS DAMPING RATE AS A FUNCTION OF GAS PRESSURE 


The damping rate of the particle motion due to the background gas is given by /q = 77 //to, where to is the particle 
mass and 77 / = is the coefficient of friction for a spherical particle of radius Vd in a fluid of dynamic viscosity 

/i. This expression is for a particle much larger than the mean free path of the gas bath Amfp- For a nanoparticle in 
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a rarefied gas, as assumed in Eq. (9) in the main text, the damping constant becomes |S6| 

^ dnuTd . hTTurd 0.619 /. O.SlKn 

In = - X correction = - 1 H- 

m m 0.619+ KnV 0.785 + 1.152Kn + Kn 



(S28) 


where the correction term is in terms of the Knudsen number Kn = Anifp/fd oc 1/ Pr^. Using the fact that Amfp ~ 70 nm 
at atmospheric pressure, and matching the experimentally measured rates |S31 IS5| . we can write the damping for a 
general pressure and particle size as 


rd 271x10® Hz/ 0.31Kn 

70 nm 0.619+ Kn V 0.785 + 1.152Kn + Kn^ 


(S29) 


In order to ensure that this model of damping holds in the deep quantum regime [S7], we consider the effects of 
individual scattering events as described by the master equation |S8| 


p = Y^VpV^ -\[v^v,p], 

V 


where the Lindblad operators are given by 


V 



(S30) 


(S31) 


where Ap is the momentum transfered to the nanoparticle during a collision with a gas molecule. Now if the momentum 
kick is sufficiently small, then we can approximate 


exp (y—iApq/h)j « 1 — iApq/h, 

which ultimately allows us to recover the Brownian contribution, derived in 


(S32) 

to our master equation presented in 


the main Letter. In order for the approximation of Eq. (S32) to hold, it is necessary that Z\p/coh/^ < 1, where £coh 
is the coherence length of the system. Considering the coherence length of the system in its ground state, £coh ~ £z- 
Using the parameters listed in the main article for the system at T = 4K, and using the maximum RMS thermal 
momentum change Ap = assuming a gas made of N 2 gives 


Apt:, 

h 


0.6 < 1 , 


(S33) 


which ensures that our linearized treatment of Brownian motion is valid over all considered regimes. This inequality 
actually overestimates the magnitude of a typical momentum kick as Hydrogen and Helium comprise the majority of 
the species of gas molecules in an environment at 4K, and are much lighter than N 2 . 


CONVERTING A MASTER EQUATION TO A LANGEVIN EQUATION 


Consider a master equation in standard Linblad form 

3 


(S34) 


where iLg is a Hamiltonian representing the unitary evolution and 'D[Lj]p is the standard Lindblad dissipation su¬ 
peroperator for the Lindblad operators Lj [S9]. Now following the quantum-state diffusion model presented by Gisin 
and Percival the density matrix can be considered as the expectation over pure states e.g. p = E [|'!/)('0|], 

whose dynamics are stochastic in nature. The (stochastic) differential increment for the equation of motion for Itjj) 
representing the unraveling of the master equation given by the QSD model is |S10I IS11| 


IdT/) = ^Hdt\^)+Y,\Y (2 {l]} L, - l]l, - {l]) {L,)) dt + {L, - {L,)) V^dlU, 




(S35) 


6 

















where dWj are complex differentials representing Wiener processes defined by E [dWj\ = E 
0, and dW* dWk = Sjk dt. Note, this complex Wiener process is equivalent to dW = (dRi 


dw: 


= E [dWj dWk] = 


idi?2)/'\/2, where dRj 

are real valued Wiener increments. Now such a model is not unique, but represents an open quantum system which 
is continuously monitored by it’s environment and thus provides an appropriate trajectory or unraveling for 

a system that we will ultimately wish to observe in terms of a PSD or measurement spectrum (i.e. which we will 
continuously monitor). 

Now we are interested in expressing the random forces acting on our system due to interaction with the environment 
(as expressed by the Lindblad operators Lj). Therefore computing the increment for the momentum p of the system 
gives 


dp = (dV'lpIV') + (V’bl db) + {dip\p\ dip) 


3 

= [P] di + X! + P ) 


pdt 


0{dt^) 

^(dB]p- 


■pdSj 


(S36) 


where dBj = {Lj — {Lj))y^dWj. Applying this to our master equation gives us our Langevin equation for the 
stochastic forces listed in Eq. (15) in the main article. 
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